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We develop an analytical tool to extract bulk electronic properties of unconventional superconductors through 
inelastic neutron scattering (INS) spectra. Since the spin excitation spectrum in the superconducting (SC) 
state originates from Bogoliubov quasiparticle scattering associated with Fermi surface nesting, its energy- 
momentum relation-the so called 'hour-glass' feature-can be inverted to reveal the Fermi momentum disper- 
sion of the single -particle spectrum as well as the corresponding SC gap function. The inversion procedure 
is analogous to the quasiparticle interference (QPI) effect in scanning tunneling microscopy (STM). Whereas 
angle-resolved photoemission spectroscopy (ARPES) and STM provide surface sensitive information, our in- 
version procedure provides bulk electronic properties. The technique is essentially model independent and can 
be applied to a wide variety of materials. 



The study of cuprate superconductivity has led to the dis- 
covery of a number of unexpected relationships between 
seemingly different spectroscopies. STM provides an out- 
standing example of this: while basically a real-space probe, 
it can extract momentum-(/c)-space information of the elec- 
tronic Fermi surface and SC pairing symmetry usually asso- 
ciated with ARPES. Elastic scattering of Cooper pairs leads 
to the QPI pattern measured by STM.|ni,f2l| By analyzing the 
QPI pattern as a function of scattering angle (in g'- space) it 
is possible to reconstruct both band structure and gap infor- 
mation in /c-space. While the technique was originally de- 
veloped for d-wave cuprates, it is finding wide applications 
in other materials, including pnictides, where the SC gap is 
probably of s± symmetry, and in topological insulators in the 
non- superconducting state. 

The development of a similar technique for INS would have 
obvious advantages. First, ARPES and STM are surface sen- 
sitive probes, which in practice means that they can only be 
performed on materials that are readily cleaved. For the same 
reason, they are mainly restricted to quasi-two-dimensional 
materials, and the question always remains how sensitive the 
results are to surface effects, such as pinning, reconstruction, 
and excess scattering. In contrast, INS is a bulk probe which 
does not need special surface preparation and can readily be 
applied to three-dimensional materials. For instance, INS re- 
sults on heavy fermion materials were available years before 
the first high-resolution ARPES studies. Here we demonstrate 
that inelastic scattering between the particle and hole Bogoli- 
ubov quasiparticles generates a similar magnetic quasiparticle 
scattering (MQPS) profile which can be probed directly by 
INS measurements. ISHIl Taking all three spectroscopies to- 
gether, we are able to establish a definite consistency between 
the r-space (STM), fc-space (ARPES), and (7-space (INS) 
dynamics of the unconventional Bogoliubov quasiparticles. 

The overall phenomenology of neutron scattering in 
cuprates is well-known, and a number of universal features 
have been experimentally identified |[6HT4l and theoretically 
interpreted.|3-5| The results indicate that a distinct low en- 
ergy magnetic mode is present near the antiferromagnetic 
nesting vector Q = (tt, tt) in almost all the cuprates. The 
intensity of this mode is enhanced in the SC state while its 



energy scales uJresiQ) 2Asc for all cuprates fT5l[T6ll , sug- 
gesting a close connection of these modes with SC pairing. 
The dispersion of spectral weight away from this resonance 
peak also has a universal character, forming an 'hour-glass' - 
like pattern centered on the resonance mode and displaying 
a 45^ rotation on passing through the resonance peak. Be- 
low the resonance energy, spectral weight disperses along the 
Cu-0 bond direction, while above the resonance the disper- 
sion peak lies along the diagonal direction. Despite an over- 
all universality, many features of this 'hourglass' dispersion 
are highly material specific and we show that the difference 
comes mainly from the nature of the pseudogap order among 
other band structure related properties. 

We begin with describing the inversion procedure in INS, 
which is analogous to the common inversion procedure per- 
formed using the QPI features in the STM spectra. Fig. 1(a) 
depicts a typical form of the Fermi surface (green shadings) 
which is assumed to be truncated below the antiferromag- 
netic Brillouin zone due to the presence of underlying pseudo- 
gap ordering. When the superconducting gap has a (i— wave 
form Afc = Ao{cos{kxa) — cos{kya))/2, where a the in- 
plane lattice constant, the spectral function at any particular 
energy u < 2Ao has a maximal intensity at eight bright- 
points K{uj) that develop on four banana-pockets which sat- 
isfy uj = Ax{oo) (red dots). In STM experiments, elastic scat- 
tering at energy uj is dominated by scattering between these 
points, leading to a form of Friedel oscillations known as QPI, 
satisfying the condition ||T]|2l 

^ ^^(a;) ^K{uj)+qi- (1) 

The long dashed cyan arrow in Fig. 1(a) illustrates one possi- 
ble -vector. 

The INS experiments, which measure the imaginary part of 
the transverse spin-spin correlation function x" ^ can also be 
understood through an analogy with the QPI pattern. 1 3 , 5 1 In- 
deed, since x is related to the joint density-of-states (JDOS), 
the neutron scattering will be dominated by transitions from 
banana tips below the Fermi level to banana tips above 
the Fermi level, as illustrated by the dotted cyan arrow in 
Fig. 1(a). The energy of a particular transition will be given 
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FIG. 1 . Origin of MQPS vectors, (a) Sketch of cuprate Fermi sur- 
face, showing four 'banana' -shaped regions (green shapes) with in- 
tense spectral features at their tips (red dots). The gap as a function 
of angle is shown on one banana. Above the Fermi pocket, the super- 
conducting gap is cut off by competing-order pseudogap (black solid 
line). On this gap curve the contrasting behavior of QPI and MQPS 
vectors is illustrated: the QPI vectors connect two banana tips at the 
same energy, as shown by the long dashed cyan line, while the MQPS 
vectors actually connect banana tips across the Fermi surface, from 
an electron tip at a; < to a hole tip at cj > 0, shown by the dotted 
cyan arrow for the Qs branch. For convenience the various MQPS 
vectors are defined by projected horizontal arrows, labelled Qi — Qs, 
illustrating the two bright points that they connect, (b) The ^- space 
representation of the MQPS. Symbols represent the Q^— vectors of 
same color as shown in (a). The associated arrows point in the direc- 
tion that the Qi vectors shift with increasing energy along the Fermi 
surface, (c) The energy and momentum dispersion of two Qi vectors 
are shown along the diagonal direction [Qx = Qy]. [The intensity 
of Qs is expected to be low (see text).] The black solid line de- 
picts the spin- waves which are gapped below the SC gap while the 
black dashed line stands for the gapless spin-wave (schematic) dis- 
persion expected in the normal state. The black dashed line is the lin- 
ear spin-wave expected in the non- superconducting region, (d) Both 
spin-wave (solid black) and MQPS (solid blue) have monotonic in- 
tensity variation as a function of energy, and as they meet at Q, the 
total intensity (red line) attains a sharp peak. 

by 

uj{Qi) = \Ak\ + |Afc+Qj = Ao [\gk\ + l^fc+gj] , (2) 

where the d— wave structure factor is = (cos (kxa) — 
cos (kya))/2. 

While the same banana points contribute to both QPI and 
neutron scattering, there are actually more gr— vectors in the 



latter case owing to more possibilities of inelastic scattering. 
In elastic scattering there are 7 g^^— vectors at any energy that 
connect a given banana point to the other banana points. dm 
However, in INS there is a coherence factor which allows only 
the poles for which and A^+q . have opposite signs, since 
the magnetic neutron scattering cross-section is odd under 
time-reversal symmetry. fT-'Sl Hence we identify four neutron 
scattering vectors which dominate in cuprates and give rise to 
a spin excitation profile, which can be called MQPS pattern, in 
analogy with the QPI patterns. The MQPS vectors can be de- 
noted as Qi^2,3 ^3,6,7. see Fig. 1(a), where the lower-case 
grs are the corresponding QPI vectors. d 12 

There is a second, more important difference from QPI. 
Note that for fixed uj Eq. [T] is the equation of a point in q- 
space, while Eq. |2]is the equation of a curve. Because there 
will be a different Qi for each pair A^ and A^+q. which 
satisfy Eq. [2] However, we show below that this is not a 
problem: the gap and FS can be reconstructed from INS data 
along any Q — uj cut. INS data are mostly available along 
the diagonal cut Qx = Qy, and along the bond direction, i.e. 
[(0,7r) (27r,0)] and equivalent cuts as a function of energy. 
We therefore define the spanning vectors as Qi and Qs along 
the diagonal direction and Q[ and along the bond direc- 
tion [note that Q2 does not lie on these two special cuts]. The 
diagonal cut is a particularly simple choice, since along this 
cut the MQPS vectors are exclusively associated with special 
points, for which A^ = — A^+q . . For these special points the 
analogy with QPI becomes exact, with ujmqps = '^^qpi- 
Along the bond direction, the situation is similar as discussed 
below. 

The Q-space positions of the poles of Eq. 2, or the MQPS 
pattern, are plotted schematically at a constant energy in 
Fig. 1(b). The associated arrows indicate the direction each 
vector moves as the excitation energy increases which is sim- 
ply a manifestation of the shape of the Fermi pocket and the 
(i— wave superconducting gap values at each Fermi momen- 
tum. To see this, we draw an illustrative energy versus mo- 
mentum dispersion relation for two Qi vectors along the diag- 
onal direction in Fig. 1(c). Branch Qi disperses from uj = 0, 
where the incommensurate Qi -vector connects pairs of nodal 
points, and hence measures the width of the hole pocket, to the 
resonance peak at (tt, tt) at the maximum uj = 2ujres connect- 
ing the hot-spots at the boundary of the magnetic Brillouin 
zone - that is, Qi is the dispersion branch associated with 
the magnetic resonance peak. Branch Qa represents an in- 
trapocket scattering, and the corresponding branch starts from 
Qs = 0, UJ = at the nodal point, dispersing towards a 
Q 3 -vector which spans the length of the hole pocket when 
UJ = 2ujres- Thus both curves attain the maximal uj when the 
banana points correspond to the hot spots along the Brillouin 
zone diagonals - the maximal diagonals of the hole pock- 
ets. A similar phenomenon is found in the experimental QPI 
spectra. lUO In summary, the inelastic MQPS intensity pat- 
tern at an energy uj = 2 A in neutron scattering exactly cor- 
responds to the elastic QPI pattern seen at an energy u = A 
in STM, and hence can also be used to reconstruct the Fermi 
surface and SC gap properties. 

To illustrate how the above results play out in practice. 
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we compare them to experiment and to realistic calcula- 
tions of INS. We calculate the INS spectrum of the cuprates 
from a one-band Hubbard model, with self-energy correc- 
tions from a GW calculation called the quasiparticle-GW 
(QP-GW) modeL|17, 18] The QP-GW hamiltonian consists 
of four components which are calculated self-consistently: 
Hlda + HsDW + Hsc + The electronic dispersion, 
Hlda, is based on a tight-binding fit to the material- specific 
first-principles single band dispersion of the antibonding com- 
bination of Cu dx2_y2 and O p^/y orbitals.| 19] The pseu- 
dogap state of the cuprates is treated as a spin-density wave 
(SDW), HsDW, based on a Hubbard term in the Hamiltonian, 
which is treated using standard random-phase- approximation 
(RPA) theory. 1 20 1 Below T^, c^-wave SC order develops Hsc 
which couples naturally to the SD W state.] 21] The result- 
ing energy spectrum is E{k) = Y^^(E^^)^(^yHKA|, where 
the non- superconducting dispersion for both spin states are 
Es = ^+{k) ± Eo{k) [for E^ = ^-{kf + G\ ^^{k) = 
(^(/c) + ^(/c + Q))/2, ^(/c) is the non-interacting band in the 
Bloch state, and G is the effective SDW gap which produces 
the pseudogap above the antiferromagnetic zone boundary]. 
The SDW gap causes a substantially reconstructed Fermi sur- 
face (FS) at low-temperature as shown in Fig. 2(c). At low 
temperature the wave superconducting gap coexists with 
the SDW state. By fixing U and the pairing interaction V, 
to account for the experimental values of pseudogap and su- 
perconducting gap respectively (see Table I), there are no free 
parameters in the susceptibility calculation. Finally, we calcu- 
late the self energy due to the magnetic and charge excitations 
which renormalizes the overall dispersion by a momentum in- 
dependent mass renormalization. 

The full susceptibility is computed in all spin channels 
within random phase approximation (RPA). In RPA, the res- 
onance pole is determined by the real part of the Lindhard 
susceptibility, Xo which attains a logarithmic divergence at all 
the Bogoliubov quasiparticle scattering vectors. Simultane- 
ously Xo possess discontinuous jump due to Kramer's Kronig 
relation. In this spirit, the positive divergence in Xo or the dis- 
continuous jump in Xo can be used as an indicator of the ob- 
served the magnetic spectra in the SC state. In the SDW state, 
the situation is more complicated, nevertheless, the overall 
phenomenology can be discussed approximately by the imag- 
inary part of the transverse spin susceptibility which mainly 
consists of three important factors (see appendix for more de- 
tails): 



Xo(g,cj) ^ ^'^'^A{k,q)C{k,q)d{uj -ujres{k,q)i3) 

k 



Here, ujres{kj q) = E{k) + E{k + g), which converges to 
Eq. 2 on the Fermi surface [E^{k) = 0]. Since the INS spec- 
trum is proportional to the RPA susceptibility x'^ the asso- 
ciated intensity of a given pole in the spectrum is controlled 
mainly by the two coherence factors A and C associated with 
antiferromagnetic zone folding and SC gap symmetry break- 



ing, which can be approximated for discussion as 

^^^'^^ - 2 E{k)E{k + q) ) ■ 

At the normal state Fermi surface [E^{k) = 0], the su- 
perconducting coherence factor reduces toC = 1/2(1 — 
sgn(^fc)sgn(^fc+q)) which attains its maximum value of 1 
whenever and A^+g have opposite signs, thereby explain- 
ing why the spectrum is dominated by the four magnetic scat- 
tering channels Qi with {i = 1 — 4). Furthermore, at small 
\q — Q\, the SDW coherence factor, Eq. 5, simplifies! 22] to 

A ^ l-(r (fe)+r(fc+9))GV4E4(fc) = l-0{{Q-q)^), 
which attains its maximum value of 1 at g = Q. This explains 
why the intensity for Qi branch being closer to Q dominates 
over the other branches and the intensity gradually increases 
to its maximum value at Q. 

Fig. 2 presents our main result, showing that the dominant 
spectral features in neutron scattering correspond exactly to 
the MPQS features, in a SDW superconductor. For a spe- 
cific example we analyze YBa2Cu306+y (YBCO). The color 
plots in Fig. 2(a) and 2(b) show two-dimensional g — in- 
tensity maps of the calculated imaginary part of the RPA 
spin susceptibility (x'O along the diagonal and the bond di- 
rections for YBCO Sit y = 0.85. Shown superimposed on 
the left sides of Fig. 2(a) are the calculated MQPS disper- 
sions, with Qi = (tt ± J, TT ± J) plotted as a dashed blue 
line and Qs = as a dashed cyan line, computed 

from Eq. 2. Similarly, in Fig. 2(b) the dashed red [yellow] 
line represents the equivalent branch along the bond direc- 
tion: Q[ = (^,^± J)/(^± J,^) [Q^ = {7r,±6)/{±6,7r)]. 
The susceptibility maps present complex patterns with numer- 
ous dispersing features. Nevertheless the MQPS branches Qi 
and Qs capture the essential shape of the two downward dis- 
persions, while a spin wave branch is clearly seen dispersing 
upward away from the magnetic resonance feature. 

There are two main sources of intensity variation for each 
Q— vector. As discussed earlier, the SDW coherence factor 
is strongly momentum dependent and reaches maximum at 
Q. Since branch Qs stops dispersing well before it reaches 
Q, its intensity is relatively low while Qi branch gains more 
intensity as it moves towards the resonance. The equivalent 
QiiQs) along the diagonal direction in Fig. 2(b) has twice as 
large intensity as QiiQs) in Fig. 2(a). This is due to an over- 
all degeneracy factor. The scattering vectors Q[ ^ have one 
commensurate direction and one incommensurate one while 
Qi^s are incommensurate along both x— and ?/— directions 
except at the resonance at Q. Thus Q[ 3 connect twice as 
many Fermi surface points as any other Qi,3- As a result the 
intensity of the magnetic spectrum along the bond direction 
[in Fig. 2(a)] is twice as large as that along the diagonal direc- 
tion [in Fig. 2(a)] [also see several intensity plots at relevant 
momentum cuts in supporing Fig. S4]. Therefore, in the su- 
perconducting region, the intensity profile is rotated along the 
bond-direction as shown in Fig. 3(a). We find good agreement 
with the experimental results, represented by symbols of vari- 
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FIG. 2. Reconstruction of Fermi surface from the INS spectra. 

(a)-(b) Computed x'{^^^) is plotted along the diagonal and bond 
directions respectively (see inset) in the superconducting state of 
YBCO. The experimental data for YBCOe+y are shown in various 
symbols [green stars {y — 0.85) =Ref. |8|, blue squares {y — 0.85) 
=Ref. nil, blue triangles {y = 0.95)=Ref. ||23l]. The data along 
the bond direction in (b) are obtained for y — 0.6 (red symbols) 
from Ref. |f5|. The blue and cyan dashed lines in (a) are the plot of 
Eq. 2 for Ql and Qs (red and yellow dashed lines are for correspond- 
ing equivalent vectors Q'l and Q3 along bond directions in (b). The 
brown symbols give the resonance spectrum calculated from Eq. 2 
using ARPES dispersions|25 1 for y — 0.6 as input [see (c)]. The 
experimental and theoretical lines are only plotted in one side of two 
equivalent directions for clarity. We plot x'{(l^^) log-scale to 
highlight weaker features. The black line is a guide to the eye for the 
dispersion expected for gapped spin- waves, which constitutes the up- 
ward dispersion in the resonance spectra, (c) Plot of computed Fermi 
surface of YBCO. The blue symbols are the locus of the Fermi mo- 
menta determined from the experimental neutron data shown in (a) 
solving Eq. 2 while red symbols are the same but along the bond di- 
rection shown in (b). The brown symbols are the ARPES Fermi sur- 
face from Ref. |25|. All symbols have four-fold symmetry, but are 
plotted here only along one particular quadrant for visualization. The 
blue and red arrows are the scattering vectors Qi and Q'l. (d) The 
extracted superconducting gap plotted as a function of Fermi surface 
angle [O'^ being the antinodal direction and 45^ the nodal direction] 
from the neutron data of (a) and (b), compared to the ARPES data 
(brown symbols) of Ref. [26 1 and the present theoretical curve. 

ous colors |[8l[T2l|23l EH, on the right-hand side Figs 2(a) and 
2(b). 

However, here our principal task is to demonstrate that the 
Fermi surface and superconducting gap can be reconstructed 
from the experimental data. The lower branch of the experi- 
mental data which disperses downward from (Q, ujres) is as- 
sociated with the Ql scattering vectors. Therefore, by solv- 



ing Eq. 2 at the experimental points, we determine a map of 
kp which agrees very well with the theoretical Fermi surface 
as well as with the ARPES FS. Note that we find the same 
Fermi surface by utilizing the diagonal cut Qi [blue dots in 
Fig. 2(c)] or the bond direction Qi [red dots]. For a cross- 
check, we use the experimental FS measured by ARPES EH 
[brown symbols in Fig. 2(c)] to compute the magnetic res- 
onance spectra by solving Eq. 2, brown symbols in Fig. 2(a). 
This agrees well with the INS data. Furthermore, the extracted 
values of hp versus uo can be used to determine the under- 
lying superconducting gap symmetry using Eq. 2 which as- 
sumes dx2_y2 — symmetry. This also agrees well with ARPES 
data|l26l and the present theory, see Fig. 2(d). The supercon- 
ducting gap can be extracted up to the edge of the magnetic 
Brillouin zone or the hot-spot [(j) ~ 15^], above which the 
pseudogap dominates in the spectra and the IMPQS features 
can no longer be followed. 

Above we have discussed the magnetic resonance peak and 
the lower half of the hourglass. The upper, rotated half of 
the hourglass has a different origin, related to spin waves. 
The black dashed line in Fig. 1(c) shows the spin- wave dis- 
persion in the non- superconducting state at half-filling, which 
disperses away from the antiferromagnetic wave vector Q = 
(tt, 7r)|20|. When superconductivity is turned on, low energy 
states near the Fermi surface are gapped (black solid line) and 
the paramagnon features at energies uj < 2 Aq are washed out. 
The black dashed lines in Figs. 2(a), 2(b) show that the present 
calculations capture this gapped spin wave spectrum. 

Combining the IMQPS and spin wave spectra. Fig. 3 com- 
pares the -evolution of our present theoretical IMQPS maps 
[Fig. 3(b)] with corresponding experimental data [Fig. 3(c)]. 
At low energies the spectra are dominated by the IMQPS vec- 
tors, and the figure also compares a QPI map [only available 
for Bi2212], in Fig. 3(c), with a neutron scattering map of the 
odd channel in YBCO, in Fig. 3(d). Three high-symmetry 
vectors (5i,2,3 (and two equivalent vectors along bond di- 
rection Ql 2) are visible which correspond to the ^3,6,7 vec- 
tors in the QPI map [denoted by circles of the same colors]. 
In the low-energy region, the magnetic scattering profile of 
Fig. 3(al) is squarish with finite intensity along all directions. 
However, the intensity is larger along the bond-directions than 
along the diagonal as discussed above and thus the profile is 
rotated along the bond directions, consistent with experiments 
shown in the corresponding lower panel. With increasing en- 
ergy, at the resonance [Fig. 3(a2)] the intensity piles up at 
(tt, tt) with tails dispersing towards q = (0,0) which were 
recently seen in Hg-based compounds. [27] Above the reso- 
nance. Fig. 3(a3), the magnetic spectrum is purely spin-wave 
based and the profile is rotated along the diagonal direction, 
as found experimentally. 

In the above calculations, we have assumed that the super- 
conductivity couples to the SDW order. When, the spin-wave 
spectrum of upward dispersion and the IMQPS of downward 
dispersion meet at the commensurate antiferromagnetic vec- 
tor Q = (7r/2, 7r/a), a resonance peak in the intensity occurs. 
Such a resonance peak is absent in a paramagnetic ground 
state. (SHU The SDW coherence factors play an important role 
in distributing the intensity over the entire magnetic spectrum. 
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FIG. 3. Comparison of MQPS maps with experimental QPI maps. (al)-(a3) Computed spectra of x'{^) are plotted in color code for 
three representative constant energy cuts, below the resonance, near the resonance and above the resonance. (bl)-(b3) The corresponding 
experimental data of YBCOe.gs from Ref. |ir| for the acoustic channel. The white arrows in the experimental curves were used in the 
discussion in Ref. 1 11 1. (c) A QPI map of the STM data|2 1 for Bi2212 is compared with the MQPS map of INS in (d). The MQPS spectrum 
is chosen in (d) for the optical channel to ease comparison with QPI maps as the magnetic profile is rotated here along the diagonal direction. 
We note that QPI are usually plotted in the range [— tt to tt] while the neutron scattering is plotted in the range [0 to 27r]; the dashed box in (c), 
(d) represents the range [0 to tt] to be compared. In (al), (bl) and (d), circles of various colors depict various MQPS vectors defined before 
while the circles with same colors in (c) compare QPI vectors with corresponding MQPS vectors. 



Eq. 4. Furthermore, the QPI-MQPS correspondence is ob- 
scured when the non- superconducting state becomes param- 
agnetic. While the Qi-vectors still play a significant role, 
the magnetic resonance peak is shifted to a lower energy at a 
pole of the dynamic susceptibility|5 |. In the magnetic super- 
conducting ground state, the susceptibility peaks correspond 
to the MPQS Qi{uj), with uj given by Eq. 2, as shown in 
Fig. 2A,2B. 

In summary, we have shown that while the spin- wave spec- 
trum of collective mode origin persists at all dopings both in 
electron and holed doped cuprates including at half-filling, it 
becomes gapped in the low-energy region u < 2A where 
the spectrum is dominated by the Bogoliubov scatterings of 
Cooper pairs. When, the spin- wave spectrum of upward dis- 
persion and the magnetic scattering of downward dispersion 
meet at the commensurate antiferromagnetic vector Q, a res- 
onance peak in the intensity occurs. The INS spectroscopy 
is the only probe which can be utilized to reconstruct the 
Fermi surface and superconducting properties of the actual 
bulk ground state. Our method of analysis is independent 
of any particular model and can be performed entirely from 
experimental inputs. Therefore, this inversion procedure will 
also be useful in further research to extract the bulk-Fermi sur- 
face topology and pairing symmetry information in newly dis- 
covered iron-selenide, 1 28 1 |29i pnictide, chalcogenide| 30 | and 
heavy fermion| 31 1 superconductors in which this information 



is still not settled. We also predict that the present formalism 
can be used to detect the electron pocket on the cuprate Fermi 
surface, if they are present in the bulk ground state. 
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Appendix A: Susceptibility in SDW+d-wave SC state 

Since the SDW state causes a unit cell doubling, the corre- 
lation functions (Lindhard susceptibilities) are tensors in mo- 
mentum space representation! 20 1 We define the susceptibili- 
ties as the standard linear response functions 

x'HQ,q',r) = ^(r.n^(r)ni^,(o)) (ai) 
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Material 


Doping (x) 


Apg (meV) 


U/t 


Asc (meV) 


Pairing Potential 










(Exp./Theory) 


(Theory) 


(Exp./Theory) 


y (meV) (Theory) 


Exp. (Theory) 


Exp. (Theory) 


NCCO 


0.15 


60 [Ret 32} 


4.1 


3.5 [Ref.l3J 


-26 


24 (37) [Ref. 33J 


4.5 (5) [Ref.. 32 J 


YBC06.85 


0.21 


50 [Ref. [341 


2.5 


35 [Ref. [231 


-98 


92 (160) [Ref. [231 


41 (40) [Ref.[8l 



TABLE I. Order parameters of the model. Experimental gap values and the resulting self-consistent values of the order parameters. Given 
the order parameters, theoretical values of ujres are calculated and compared with experiment. The value of U/t (where t is nearest-neighbor 
hopping parameter) is chosen to reproduce the experimental peudogap (Apg) whose values are presented here along the 'hot-spot' direction 
in the electron doped case and the antinodal direction for all hole doped cuprates. Similarly, the parameter value of pairing potential V is 
taken to reproduce the superconducting gap (A), whose maximum value along the antinodal direction is presented here. 1 35 1 Our mean-field 
calculations overestimate the values of Tc, presumably due to the neglect of phase fluctuations. 1.3 6 J 



where the response operators (H) for the charge and spin den- 
sity correlations respectively are 

Pq(^) = and 

k,a 

Si^r) = c^,^^(TX,Cfc,,(r). (A2) 

The cr* represent two dimensional Pauli matrices along i^^ 
direction. Cfc,cr(c^^) is the destruction (creation) operator 
of an electronic state at momentum k and spin cr. For 
transverse spin response = Sx ^ iSy whereas longi- 
tudinal fluctuations are along the z— direction only. In the 
present (tt, 7r)-commensurate state, charge- and longitudinal 
spin-fluctuations become coupled at finite doping. In common 
practice the transverse, longitudinal spin and charge suscepti- 
bilities are denoted as , X^^ X^^ respectively. We col- 
lect all the terms into a single notation as where a = a 
gives the charge and longitudinal components and a = —a 
stands for the transverse component. For the pure SDW state 
Eq. Al can be evaluated rigorously. Here we generalize ear- 
lier calculations! 20 1 for realistic cuprate band structures. For 
the combined SDW+d-SC state 



k,n,s 



(A3) 



1 ' ^ 

= N^ ^''^^ ^ C™,xr.'(fc,9,c^). (A4) 



We obtain Eq. A4 from Eq. A3 after performing the Matsub- 
ara summation over n. G is the 4 x 4 single-particle Green's 
function in the Nambu space, constructed from the Hamilto- 
nian given in the main text. The summation indices I'ii'') — ± 
gives two split SDW bands. Here, the coherence factor due to 
SDW order in the particle-hole channel is. 



2 V^ofc 
The SC coherence factors are 




Ofc+q 



(A5) 



a 



2/3 



c}, 



1 ± 



EY EY I 

k k+q 



(A6) 



Lastly the index m represents the summation over three polar- 
ization bubbles related to the quasiparticle scattering (m = 1), 
quasiparticle pair creation (m = 2) and pair annihilation 
(m = 3), as defined by 



xly{k,q,u) = 



f{K)-fiE^k^,) 



(A7) 



(A8) 



In the RPA model, the 2 x 2 susceptibility is obtained from 
the standard formulal.2QJ 



m—l 



XrPA,22{Q^^) 
X'rpA,12/2i(Q^^) 



[l+aaUx^Uq,^)]x^!{q,^) + U[x^,i{q,u;)y 



[l-[/x?f(g,w)] [1- 


^aaUx^i{q,uJ 


)] +aa[Ux'[i{q,uj)f 


[1 - UxfHQ, 


i^)]x2iiq,^) + 


c/[x?f(g,^)]' 


[l-Ux^!{q,UJ)] [1- 


^(TaUx^i{q,uj 


)]+(Ta[Ux^,i{q,oj)f 




Xii{q,to) 




[1-Ux^!{q,uj)] [1- 


^(TaUxM{q,uj 


)]+aa[Ux'[!{q,i^)f 



(A9) 

(AlO) 
(All) 



In the longitudinal and charge channel (cr = cr), the RPA cor- rections do not introduce any pole and thus all the normal 
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(0,0) Qx=Qy {2nla,2nla) 
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5- 



~ " X 

~~ X"rpa . 




(c) 









25 
03 (eV) 




"0 

m6 



50 Q 



(0,0) Q=Qy (2n/a,2n/a) 



FIG. 4. The real and imaginary part of bare susceptibility is shown 
in (A) and (B) respectively for YBCO. (D), Corresponding RPA re- 
sult plotted in a logarithmic intensity scale as in the case of Fig. 2 
of the main text. (E), All three susceptibilities are shown at the com- 
mensurate momenta cut 



structure lies above the charge gap in the particle-hole contin- 
uum. Along the transverse direction (a = —a), a linear spin- 
wave dispersion develops in the normal state which extends 
to zero energy at Q = (tt, ttI. ISqI The necessary condition to 
yield a gapless Goldstone mode is that Eqs. A9||A11 reduce to 
the self-consistent SDW order parameter, G dXq = Q, which 



is indeed the case in the normal state. 

In the SC state, this zero energy spin-wave shifts to cj = 
I Afc^ I + I Afc^+q|, due to t he pa rticle-particle (and hole-hole) 



scattering terms x in Eq. A8 These terms have finite inten- 



sity only if the SC gap changes sign at the 'hot-spot' gr,| 3 1 see 
Eq. |A6 Above the SC gap, the spin- wave term coming from 
is turned on. The crossover between them creats the 



A7 



Eq. 

'hor-glass' pattern presented in Fig. 2. 



Appendix B: Mechanism of MQPS 

The result shown in Fig. 2 of the main text is obtained from 
a coexisting state of SDW and (i— wave SC order within RPA 
framework. To see how one can use the behavior of imaginary 
part of the bare susceptibility as an indicator for the observed 
resonance behavior, we decompose the resonance spectra of 
Fig. 2(b) into its bare components as shown in Fig. 4. As 
shown in Fig. 4(a) and 4(b), the real and imaginary part of 
the bare susceptibility are related to each other by Kramers 
Kronig relation. Where the real part obtains a logarithmic di- 
vergence, see blue line in Fig. 4(c), the corresponding imagi- 
nary part possesses a discontinuous jump at the same location 
[green dashed line in Fig. 4(c)]. Within RPA, a resonance is 
possible when Xo = ^/U condition is satisfied. In the region 
where Xo is greater than zero and also attains a divergence, 
a resonance can occur for a large range oiU. In this spirit, 
a true resonance spectrum within RPA can simply be identi- 
fied by tracing the divergences in Xo or by tracking the sudden 
peaks in Xo • We emphasize that this argument holds even for 
multiband pnictide supercondcutors. I,3QJ 
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